Numerical study of SUPG and LPS methods combined with higher order variational time discretization schemes applied to time-dependent convection-diffusion-reaction equations

This paper considers the numerical solution of time-dependent convection-diffusion-reaction equations. We shall employ combinations of streamline-upwind Petrov-Galerkin (SUPG) and local projection stabilization (LPS) methods in space with the higher order variational time discretization schemes. In particular, we consider time discretizations by discontinuous Galerkin (dG) methods and continuous Galerkin-Petrov (cGP) methods. Several numerical tests have been performed to assess the accuracy of combinations of spatial and temporal discretization schemes. Furthermore, the dependence of the results on the stabilization parameters of the spatial discretizations are discussed. Finally the long-time behavior of overshoots and undershoots is investigated

Numerical studies of higher order variational time stepping schemes for evolutionary Navier--Stokes equations

We present in this paper numerical studies of higher order variational time stepping schemes com-bined with finite element methods for simulations of the evolutionary Navier--Stokes equations. In particular, conforming inf-sup stable pairs of finite element spaces for approximating velocity and pressure are used as spatial discretization while continuous Galerkin--Petrov methods (cGP) and discontinuous Galerkin (dG) methods are applied as higher order variational time discretizations. Numerical results for the well-known problem of incompressible flows around a circle will be presented

Higher order continuous Galerkin--Petrov time stepping schemes for transient convection-diffusion-reaction equations

We present the analysis for the higher order continuous Galerkin--Petrov (cGP) time discretization schemes in combination with the one-level local projection stabilization in space applied to time-dependent convection-diffusion-reaction problems. Optimal a-priori error estimates will be proved. Numerical studies support the theoretical results. Furthermore, a numerical comparison between continuous Galerkin--Petrov and discontinuous Galerkin time discretization schemes will be given

Numerical study of SUPG and LPS methods combined with higher order variational time discretization schemes applied to time-dependent convection-diffusion-reaction equations

This paper considers the numerical solution of time-dependent convection-diffusion-reaction equations. We shall employ combinations of streamline-upwind Petrov-Galerkin (SUPG) and local projection stabilization (LPS) methods in space with the higher order variational time discretization schemes. In particular, we consider time discretizations by discontinuous Galerkin (dG) methods and continuous Galerkin-Petrov (cGP) methods. Several numerical tests have been performed to assess the accuracy of combinations of spatial and temporal discretization schemes. Furthermore, the dependence of the results on the stabilization parameters of the spatial discretizations are discussed. Finally the long-time behavior of overshoots and undershoots is investigated

An optimal three-point eighth-order iterative method without memory for solving nonlinear equations with its dynamics

We present a three-point iterative method without memory for solving nonlinear equations in one variable. The proposed method provides convergence order eight with four function evaluations per iteration. Hence, it possesses a very high computational efficiency and supports Kung and Traub's conjecture. The construction, the convergence analysis, and the numerical implementation of the method will be presented. Using several test problems, the proposed method will be compared with existing methods of convergence order eight concerning accuracy and basin of attraction. Furthermore, some measures are used to judge methods with respect to their performance in finding the basin of attraction.Comment: arXiv admin note: substantial text overlap with arXiv:1508.0174

Non-nested multi-grid solvers for mixed divergence-free Scott-Vogelius discretizations

Studying high-dimensional Hamiltonian systems with microstructure, it is an important and challenging problem to identify reduced macroscopic models that describe some effective dynamics on large spatial and temporal scales. This paper concerns the question how reasonable macroscopic Lagrangian and Hamiltonian structures can by derived from the microscopic system. In the first part we develop a general approach to this problem by considering non-canonical Hamiltonian structures on the tangent bundle. This approach can be applied to all Hamiltonian lattices (or Hamiltonian PDEs) and involves three building blocks: (i) the embedding of the microscopic system, (ii) an invertible two-scale transformation that encodes the underlying scaling of space and time, (iii) an elementary model reduction that is based on a Principle of Consistent Expansions. In the second part we exemplify the reduction approach and derive various reduced PDE models for the atomic chain. The reduced equations are either related to long wave-length motion or describe the macroscopic modulation of an oscillatory microstructure

Robust arbitrary order mixed finite element methods for the incompressible Stokes equations

Standard mixed finite element methods for the incompressible Navier-Stokes equations that relax the divergence constraint are not robust against large irrotational forces in the momentum balance and the velocity error depends on the continuous pressure. This robustness issue can be completely cured by using divergence-free mixed finite elements which deliver pressure-independent velocity error estimates. However, the construction of H1-conforming, divergence-free mixed finite element methods is rather difficult. Instead, we present a novel approach for the construction of arbitrary order mixed finite element methods which deliver pressure-independent velocity errors. The approach does not change the trial functions but replaces discretely divergence-free test functions in some operators of the weak formulation by divergence-free ones. This modification is applied to inf-sup stable conforming and nonconforming mixed finite element methods of arbitrary order in two and three dimensions. Optimal estimates for the incompressible Stokes equations are proved for the H1 and L2 errors of the velocity and the L2 error of the pressure. Moreover, both velocity errors are pressure-independent, demonstrating the improved robustness. Several numerical examples illustrate the results

Growth of surface undulations at the Rosensweig instability

We investigate the growth of a pattern of liquid crests emerging in a layer of magnetic liquid when subjected to a magnetic field oriented normally to the fluid surface. After a steplike increase of the magnetic field, the temporal evolution of the pattern amplitude is measured by means of a Hall-sensor array. The extracted growth rate is compared with predictions from linear stability analysis by taking into account the proper nonlinear magnetization curve M(H). The remaining discrepancy can be resolved by numerical calculations via the finite-element method. By starting with a finite surface perturbation, it can reproduce the temporal evolution of the pattern amplitude and the growth rate. The investigations are performed for two magnetic liquids, one with low and one with high viscosity.Comment: 12 pages, 12 figure